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Abstract 

> 

OO ' We present a recently introduced real space renormalization group 

CN \ (RG) approach to the study of surface growth. The method permits us 

to obtain the properties of the KPZ strong coupling fixed point, which 
is not accessible to standard perturbative field theory approaches. Us- 
qq • ing this method, and with the aid of small Monte Carlo calculations 

^\ , for systems of linear size 2 and 4, we calculate the roughness exponent 

in dimensions up to d=8. The results agree with the known numerical 
values with good accuracy. Furthermore, the method permits us to 
£^ ■ predict the absence of an upper critical dimension for KPZ contrarily 

to recent claims. The RG scheme is applied to other growth models in 
different universality classes and reproduces very well all the observed 
Q , phenomenology and numerical results. Intended as a sort of finite size 

O ■ scaling method, the new scheme may simplify in some cases from a 

computational point of view the calculation of scaling exponents of 
growth processes. 
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The study of the non equilibrium dynamics of rough surfaces and in- 
terfaces is a topic of growing interest [| Many efforts have been 
devoted to single out the behaviors shared by apparently different growing 
phenomena. In this context the search of universality classes, permitting to 
categorize different models and systems of rough surfaces sharing the same 
scaling properties, is a central task. The main milestone in this direction 
is the Kardar-Parisi-Zhang equation (KPZ), which is the minimal Langevin 
equation capturing the essential physics of many different growing surfaces 
beyond the Gaussian linear (Edward- Wilkinson) theory |l], ^, |3| . It reads ||] 

dJ ^ = ^ 2 h + \{Vhf + V ( x ,t). (1) 
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where h(x,t) is a height variable at time t and position x in a d-dimensional 
substrate of linear size L. v and A are constants and r\ is a Gaussian white 
noise. The KPZ equation has become a paradigm of surface growth pro- 
cesses and an overwhelming amount of both numerical and analytical studies 
have been performed since it was first proposed. It is also related to other 
physical systems as: directed polymers in random media [||, systems with 
multiplicative noise ||, and the Burgers equation B. 

A central quantity of interest in the study of surfaces is the roughness 
W(L,t), defined as the fluctuation of the height variable around its mean 
value, h(t); this is: 

W\L 1 t) = ±- d Y d Hx,t)-h(t)¥ (2) 
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In many seemingly unrelated problems the large scale behavior of the rough- 
ness is observed to be scale invariant, in the sense that it grows in time as 
W(L, t) ~ t 13 until it saturates to a stationary value at a characteristic time, 
t s ~ L z , with universal exponents (3 and z. Once in the saturated regime, the 
roughness obeys W(L) ~ L a . These critical exponents are not independent, 
in fact, from a trivial scaling argument a = (3z, and as a consequence of a 
tilting invariance |], [I], a + z = 2. There is therefore only one independent 
exponent, say a. 

Standard field theoretical RG studies of the KPZ equation predict the 
presence of a roughening transition above d = 2 @; i.e. there are two 
RG attractive fixed points: A trivial (Gaussian) one with a = which 
describes a flat phase, and one describing the rough phase. This last one 
is non-perturbative, therefore standard perturbative methods fail to give 
any prediction for the exponents in the rough phase. One fundamental 
question, presently under debate, is to elucidate if there is a finite upper 
critical dimension ||. Summing up, any new systematic approach to the 
problem of surface growth would be certainly desirable. 

In this paper we review a recently introduced real-space, non-perturbative 



renormalization group scheme [|9|, 10, 11] permiting to determine roughness 



exponents in growth processes (even in strong coupling regimes) with the 
aid of relatively small Monte Carlo computations. This method can be seen 
as a finite size scaling method, permitting to determine from the analysis 
of very small system sizes (which is very convenient from a computational 
point of view) scaling exponents with reasonable accuracy, and constitutes 
a powerful tool in the analysis of growth processes. 
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Let us now describe the method as applied to a microscopic model in 
the KPZ universality class ||: The solid-on-solid (SOS) model (see 0, ||, ||] 
for the definition and details). Instead of studying its microscopic dynamics 
defined over cells of unit size Lq = 1, and height ho in a substrate of total 
length L, we describe the model at a coarse-grained level |J. For that 
purpose we first define blocks or cells, of linear size L k = 2 k and height h k = 
ho(2 K ) a . a is for the time being an unknown exponent to be determined 
latter in a self-consistent way. A system configuration is specified by giving 
the heights h(i) at scale k, at every substrate site, i (i = 1, where 
b = L/2 k ). Two configurations are considered equivalent if they differ by a 
constant height at every site. Having defined the blocks, now we introduce 
an effective coarse-grained dynamics, or growth rules at a generic scale. This 
is done by giving the transition rates from a generic surface configuration 
to a new one with one additional grown block. The transition rates should 
contain the relevant information of the microscopic dynamics. The growing 
rate at site i is defined as 

P [h(i) h(i)+h k ] = L k + X k J2 max[0, h(j) - h(i)]. (3) 

jnni 

The term L k is the contribution of the vertical growth (i.e. random de- 
position) and the sum over neighbor block sites j accounts for the lateral 
growth. We hypothesize that these are the two only possible types of growth 
in the dynamics of the SOS coarse-grained dynamics. The ratio x k = \ k /L k 
represents the relative importance of the lateral growth from an elementary 
lateral step in comparison with the growth through random deposition at 
scale k. The parameter x k , necessary to fully specify the dynamics is to be 
determined self-consistently by imposing scale invariance. 

Iterating in time the aforementioned growth rules, a statistically station- 
ary state should be reached, the averaged width of which could be deter- 
mined as a function of L k , x k and h k . At this point we can pose ourselves 
the following question: given a description of the system at scale k, de- 
fined by parameters (L k , h k , x k ), what will be the values of h k +i and x k+ \ 
at the scale Lfc+i? To answer this question we consider a given system of 
size L = L k+ 2 and describe it at two different levels of coarse graining: 
using b blocks of size L k +i = L k +2/b and then dividing each component 
block of size L k+ \ again in b blocks of size L k , or alternatively consider- 
ing directly b 2 blocks of size L k = L k+ i/b = L k+ 2/b 2 - The condition that 
the total roughness takes the same value when calculated using blocks of 
two different scales gives a renormalization equation in a way to be ex- 
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plained in what follows. In particular, we use the following a property of 
the roughness, W 2 : Partitioning the system in cells of size a, one obtains 
W 2 (L) = W 2 (a) + W 2 (L/a), i.e. the total roughness, W 2 (L), is equal to 
the roughness of the system measured using cells of size a, W 2 (L/a), plus 
the internal roughness of one of these cells, W 2 (a). First we apply the last 
equation to the case in which L = L k+ \ and a = L k = L k +i/b, obtaining: 
W 2 (L k+ i) = W 2 (L k ) + W 2 (b;x k ) where W 2 (b;xk) is the stationary rough- 
ness for b blocks for x k , which is proportional to h 2 ,, and therefore can be 
formally written as W 2 (b; x k ) = R{xk, b)h\. As shown in |J h\ = c* W 2 (L k ) 
where c is a constant [|l2[; this simply reflects the fact that the height h k at 
scale k has to be taken proportional to the roughness at that scale. Com- 
bining the last two results, we can write 

W 2 (L k+l ) = W 2 (L k )(l + cR(x k ,b)). (4) 

Applying twice this relation we can relate scales k + 2 and k, 

W 2 (L k+2 ) = W 2 {L k ){\ + cR(x k+1 , b))(l + cR(x k , &)). (5) 

Another possibility, as said previously, is to relate directly scales k + 2 and 
k, applying once the relation Eq.(||) replacing b by b 2 . This leads to: 

W 2 (L k+2 ) = W 2 (L k )(l + cR(x k , b 2 )). (6) 

Equating the r.h.s. of Eq.(]5|) and Eq.(^) we get the renormahzation equation: 

(1 + cR(x k , b 2 )) = (1 + cR(x k+1 , 6))(1 + cR(x k , b)) (7) 

which gives, once the function R is known, the flow of the renormalized 
parameter x k under changes of scale, and therefore permits one to evaluate 
the fixed point value x* that characterizes the scale invariant dynamics. 
Knowing x* , it is a straightforward task to determine the roughness exponent 
a: 

° = £?o \ pfe=r ) = \ l0fc(1 + cR{ *"- *»■ (8) 

The previously exposed method was applied in [|| to study KPZ growth 
in finite dimensions (from d = 1 to d = 8) by analyzing the case 6 = 2. The 
saturated width (i.e. the function R) was evaluated through Monte Carlo 
simulations of the effective dynamics in these small cells (namely 2 d and 4 d ) 
for different values of x. In all dimensions a stable fixed point is found with 
the associated values of a reported in table (1). 
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d 


1 


2 


3 


4 


5 


6 


7 


8 


org 

O^num 


0.502 
0.5 


0.360 
0.387 


0.284 
0.305 


0.238 
0.261 


0.205 
0.193 


0.182 
0.18 


0.162 
0.15 


0.150 



Table 1: Results of the renormalization group calculation (first row) com- 
pared with numerical results (second row). 



Observe the good agreement between numerical results obtained after 



extensive simulations [13], and our estimations using the new RG scheme 



with small cells H. These results support both the validity of the numerical 



simulations [13] and of our RG method. In particular, for dimensions larger 
than d = 4 exponents differing from the mean field prediction a = are 
observed, providing evidence that d = 4 is not the upper critical dimension 
of KPZ. 

The method can be applied to study analytically the infinite dimension 
limit |pX| 1. For that, one just has to evaluate the width by using some 
theoretical argument. Assuming that the lateral growth rate is very large 
for large dimensions (assumption that can be verified a posteriori), it is a 
good approximation (as discussed in |T(| ) to consider only surfaces with two 
layers: high and low respectively. A state of the system can be characterized 
by the number of sites in high (low) position, and a master equation giving 
the flow between different system states can be written down. Calculating 
analytically the width of its associated stationary state with the previous 
hypothesis for two different cell sizes, it is a matter of simple algebra to 
obtain a lower bound for a for large d [10] 

and the associated fixed point is stable for all d. This result excludes the 
possibility of having a finite upper critical dimension for KPZ-like processes. 

Even though we have presented the method as a scheme to analyze KPZ- 
like processes the underlying ideas are quite general and can in principle be 
applied to other growth processes. In particular we have studied linear 
growth models (i.e. KPZ with A = 0), that should belong in the free (Gaus- 
sian) universality class, and verified that, as expected |J, the non-trivial 
fixed point seems to disappear for d > 2. We have also applied the method 
to other highly non-trivial growth models that exhibit roughening transi- 
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Figure 1: RG flow for the two parametric model described in the text. 



tions even in one dimension [14]. In particular we have constructed a two- 
parameter representation of the model presented in [14]: with probability 
p a particle is added at a randomly chosen site, while with complementary 
probability X — p the value of the height at that site is changed by averaging 
it with its nearest neighbors weighted with a parameter a [11]- This model 
exhibits a rich phenomenology as the parameters, a and p are varied ||TTf| . 
There are three different scaling regimes: (i) a = X for a = X, (ii) for a < 1 
and small values of p there is a smooth phase with a = 0, and (iii) for a < 1 
and larger values of p, a = 0.5. This complex behavior is perfectly repro- 
duced when the model is analyzed the RG scheme. In particular in figure 1, 
we show the RG flow lines in the two-dimensional parameter space showing 
the presence of different fixed points with the expected values of a, and a 
separatrix signaling the roughening transition. 

All the previous results support the validity of the new non-perturbative 
real space RG scheme as an useful tool to analyze growing surfaces. The 
method itself can be considered as a sort of finite-size scaling procedure, 
and provides therefore a computational tool to facilitate the computation of 
scaling exponents. 

We acknowledge interesting discussions with A. Maritan, G. Parisi, A. 
Stella, C. Tebaldi and A. Vespignani. This work has been partially sup- 
ported by the European network contract FMRXCT980183, and by a M. 
Curie fellowship, ERBFMBICT960925, to M.A.M. 
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